One loop matching coefficients for a variant overlap action— and some of its simpler 
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I present one-loop perturbative calculations of matching coefficients between matrix elements 
in continuum regulated QCD and lattice QCD with overlap fermions, with emphasis a recently- 
proposed variant discretization of the overlap. These fermions have extended ("fat link") gauge 
connections. The scale for evaluation of the running coupling constant (in the context of the Lepage- 
Mackenzie fixing scheme) is also given. A variety of results (for additive mass renormalization, 
local currents, and some non-penguin four-fermion operators) for naive, Wilson, clover, and overlap 
actions are shown. 



O . I. INTRODUCTION 

o 

This paper presents one-loop perturbative calculations of matching coefficients between matrix elements computed 
using a spatial lattice regulator (and measured in lattice simulations) and their equivalent continuum-regulated values. 
The principal lattice fermion action studied is a recent formulation of overlap fermions [1]) built from a "kernel action" 
with nearest and next-nearest neighbor fermionic interactions and "fat link" gauge connections [2], specifically HYP- 
blocked [3] links. Along the way, I present some new results for a number of simpler actions-naive, Wilson, and clover 
fermions with fat links, (including a few results for APE-blocked [4] links), and overlap actions with Wilson or clover 
action kernels. I also present results for the scale at which the running coupling constant is evaluated (the so-called 
f~>) ■ q* value), using the scheme of Lepage and Mackenzie [5]. 

| Fat link actions were originally developed [6] in order to improve the chiral properties of lattice fermions. In 
£\J . the context of perturbation theory, "chiral improvement" for nonchiral actions means having small additive mass 
renormalization, vector and axial current renormalization factors nearly equal, scalar and pseudoscalar factors also 
nearly equal, and suppressed mixing of four-fermion operators between different chiral sectors. Fat link actions do 
all these things. I have found that using these actions as kernels for overlap actions seems to have computational 
advantages over using thin link actions. Overlap actions with fat links have matching factors which are much closer 
to unity than overlap actions built of thin link fermion action kernels. 

Several authors have published calculations of lattice-to-continuum matching factors for a variety of processes for 
the overlap action with the thin link Wilson action as its kernel (for a partial list, see [7-11]). Perturbative calculations 
for more complicated actions is in principle not much more difficult to do, once the Feynman rules are constructed. 
One just has to organize the calculation, and not carry it too far. Intermediate results for matching coefficients 
^ ' (particularly Table II) may be a more important part of the paper than the tables of results. These formulas are 
certainly not new, but in most published papers they appear embedded in action-specific discussions, and they are 
often hard to find until after one has rederived them. 

Whether or not perturbation theory can - or should be - used to convert lattice numbers to continuum ones is a 
matter of debate. One sometimes sees blanket statements about the use of perturbation theory in lattice calculations, 
as if perturbative calculations per se were unaffected by the particular choice of lattice action. But clearly different 
actions have different properties when the cutoff is not taken away. At nonzero cutoff, it is a practical question for 
any particular action, to ask how well perturbative calculations perform. At sufficiently small renormalized coupling 
the matching from one scheme to another can be done with small uncertainty, the uncertainty being due to higher 
order terms in a perturbative expansion in the renormalized coupling. 

To be specific, let's write a one loop matching coefficient for an operator which does not undergo mixing as 
Z = 1 + a s /(47r)A In today's simulations, typical values of a s /(47r) ~ 0.01 — 0.015. With standard actions, one often 
finds that the A's are large. For standard thin-link actions A is order 10-40, so the matching factor differs from unity 
by an amount of 0.1 to 0.6. In the actions considered here, typically the A's are smaller than about 5 in magnitude, 
corresponding to a shift from unity of 0.05-0.07. One presumes, therefore, that the perturbative expansion of Z is 
better behaved. One's results are also less sensitive to the choice of scale q for a s (q 2 ) if the coefficient is small. 

A brief report of perturbation theory for Wilson and clover actions with APE-blocked links has been presented in 
Ref. [12]. The qualitative features of the work done here are all anticipated in that paper. Here I will focus mainly on 
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HYP blocking simply because it gives more smoothing without delocalizing the action. Applications to the overlap 
in this work arc also new. Results for staggered fermions using HYP links have recently been presented by Sharpe 
and Lee [13], and Lee [14] has discussed general features of perturbation theory for fat links. 

In Sec. II, I describe the simple ingredients I used for lattice perturbation theory. Results of selected calculations 
for currents are presented in Sec. III. Appendices contain Feynman rules for the actions studied. 



II. INGREDIENTS 



A. Matching and Scale Setting 

A "typical" matching coefficient between a lattice regulated quantity and a dimensionally-regulated quantity (mod- 
ified minimal subtraction, MS, for example), for an operator which does not undergo mixing, is Z = 1 + (g 2 / (\(>it 2 ))Z 
where 

Z = I JTs^ Ilatt - W 

l£— is the finite part of the continuum 2lu = 4 — 2e dimensional integral and I latt is a lattice integral. For a process 
with Nf external fcrmion legs, /, = Tj — Ti\Nf/2 where Tj is the vertex renormalization and Si is the fermion wave 
function renormalization. All these quantities will be evaluated at one loop. In this work the internal momentum 
integration variable will always label the momentum flowing through the gluon line. I will always work at zero 
momentum for external particles and fcrmion mass m — 0, regulating any IR (infared) divergence of the diagram with 
a gluon mass A. Of course, there is nothing deep about these choices; they are made purely for expediency. I-^ s will 
take the generic form 

Ims = / ^^f \ 2[ J +x2) (A + Be) = A{\ - lE+ log( 47 r)} + Alog ^ + A + B. (2) 



IE. 

variables by appropriate powers of the lattice spacing and write 



The term in curly brackets is simply discarded to give I-jjg- In the lattice integral, we can scale all dimensionful 



Ilatt = 16tt 2 / I(ak, ap, am, aX) (3) 

Jak 

where J ak = Ylj Jl^ d(akj)/(2ir) will be the symbol for integration over the (rescaled) momentum hypercube. 

If Ijjg has an yllog(^ 2 /A 2 ) term, Ii aU will have an Alog(l/(A 2 a 2 )) IR divergence, too. It can be separated out by 
writing the integrand as 

„ 2 f, T „ „ „ 6»(7r 2 -fc 2 ) x 2 f d 4 k t 6(-K 2 -k 2 ) 

I latt = l6n jmap,am,aX)-A k2{k2 + a2x2) ) + 16, J ^ fc2(fc2 + ^2) 



= J + Aloe 



TT 2 

a 2 A 2 



(4) 



The first term of Eq. 4 is IR finite, and one can set A = in it. Thus 

Z = Alog(/i 2 a 2 ) + A(l - log, 2 ) + B-J. (5) 

This would be the end of the story if we did not want to choose a scale q* for the coupling constant in Eq. 1. This is 
(unfortunately) a problem whose solution involves at least a two-loop calculation (see the discussion in the Appendix 
of Ref. [15]). Absent such a calculation, one could either let q* vary over some "reasonable" range (1 < q*a < tt at 
lattice spacing a, for example), or construct some physically motivated ansatz for the scale. The best known example 
of such a construction is that of Lepage and Mackenzie [5] : Imagine that one has some process parameterized by a 
one loop integral, and assume that its higher order behavior is dominated by gluonic vacuum polarization, so 



I(q)d /i q -> a{q*) / I^q w / a{q)I(q)d\. (6) 
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Expanding the coupling a(q) — a(q*) — f3 a(q*) 2 log(q 2 / q* 2 ) + . . ., self consistency requires that the coefficient of 
a(q*) 2 vanish, or 

«'->-' £ T^^- !-«««■>» < 7 » 

For the actions studied here, it often happens that the Z coefficient of Eq. (f ) is close to zero, and the calculation of 
q* using Eq. (7) produces absurd results. In that case, I substitute the higher order expression of Hornbostel, Lepage 
and Morningstar [16], 

log(g* 2 ) = ((log( g 2 )))±[-a 2 ] 1 /2 (8 ) 

with 

a 2 = ((logV)» - «log(? 2 ))) 2 (9) 

and 

« logV )» - S -^P - 1 

is the weighted average analogous to Eq. (7). 

This is all well-defined for finite pure lattice expressions, but when the operator has an anomalous dimension, it 
is not obvious what to do. In that case, I evaluate Eq. 7 following a prescription learned from C. Bernard [17]. 
His proposal is to construct a combination of dimensionally-regulated integrals whose sum gives the term in curly 
brackets, and to subtract them from the integral Eq. (2) to produce a finite integral in 4 — 2e dimensions. The e — > 
limit of the subtracted integral can then be taken, leaving a UV-finite four-dimensional integral for r^-^ which can 
be combined with the integrand of Eq. 4. To do this, consider the two integrals 

/ij 2u> k 1 1 

J^(P 2 r k2{k2+ ^ (A + Be) = A{-- lE+ log( 47 r)} + A + B (11) 

and 

'2 = 16^/ ^(^ W ^(A + Be)=A{^- lE + lo g (^)} + B. (12) 

If a + b = 1 and a(A + B) + bB = 0, al\ + bl 2 is equal to the difference T-p-j — ^jjg an d we can combine the three 
expressions under one integral, then take the e — > limit, to write 

I^g = AJi + BJ 2 (13) 

where, pushing the IR divergence from the MS integral into the lattice integral 

2 f d 4 k r l-6>(7r 2 -fc 2 ) 1 1 

Jl ~ bn J J2^ [ k 2 (k 2 + X 2 ) " 

and 

' h = 16n2 Jw? [ k 2 (k 2 + v 2 ) WT7Y ] - (15) 

The integrands of Ji, J 2 , and J are then used in Eqs. 7 or 8 (with appropriate powers of log(g 2 ), of course). Naturally, 
the particular choice of integrals J\ and J 2 are not unique, but because these integrals are typically small in magnitude 
this is not generally a practical problem. Notice that the q* scale depends on fj,. I will present results only for the 
case fia=l. 

Presumably other prescriptions can be devised. Their results will probably only depend on the coefficients A and 
B in Eq. 2, so I will tabulate those parameters below for the processes considered. 

The definition of q* for matching coefficients for operators which mix, like the electroweak penguin operators, can 
be made as follows:. Eq. 7 expands into a matrix equation 
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\ogq* 2 Z = £ 1 



(16) 



where now Z and L\ are matrices. There is a basis in which log q* 2 is diagonal, found by solving the eigenvalue equation 
CiZ^ 1 — 1 logg* 2 = 0, and in the original operator basis, the matrix of scales is logg* 2 = Wilog q* 2 )di ag W~ 1 , where 
W is the basis-transformation matrix. Each eigenvector's scale can be translated into a coupling, and 

Zij = Sij + W[a s (q*)/(4Tr)] diag W- 1 Z. (17) 

The higher order formula Eq. 8 does not have an obvious matrix transformation. One possibility is to rotate the 
higher order matrix £ 2 into the basis which diagonalizes L\Z~ X , and use the diagonal entries in Eq. 8. 

This is unsatisfactory, but again, the small size of Z and the d 's for fat link actions means that the spread in the 
Z— factor remains small as the coupling is varied. 



B. Actions and Feynman rules 

Feynman rules for ordinary (non-overlap) discretizations of the Dirac operator can be constructed using standard 
techniques. I will parameterize the free massive Dirac operator as 

d(p, to) = hf,Pn(p) + A(p) + to; (18) 

I implicitly assume the normalization that at small p, p^(p) — p^ and X(p) = 0(p 2 ). (This factorization might not 
be appropriate for an approximate fixed-point action, designed to follow out a renormalized trajectory as the mass is 
varied: p M and A would both be functions of the bare mass [19].) Of course, the propagator is the inverse of d. 

While I will not specify an explicit form for the gluon propagator in any expression, I will only present results for 
the Wilson gauge action. Its propagator includes an IR regulator mass A and gauge parameter £ 

Giaf ¥Tx~ 2 ( 9) 

with = 2/asin(fc AJ a/2) and k 2 = ^ k 2 . Propagators for other actions can be constructed by inverting the gauge 
field's linearized equation of motion. It is convenient to be able to vary the choice of gauge to test results. 

The fermion and gluon propagators are of course diagonal in color space. 

The massless overlap operator is defined so that its eigenvalues lie on a circle of radius x , so 

D(0) = x {l + (20) 



where z = d{—x^)/xQ = (d— xq)/xq and d(m) = d + in is the massive Dirac operator for mass to (i.e. xq is equivalent 
to a negative mass term and d = i'f^p^ + A as above.) The overall multiplicative factor of xq is a useful convention; 
when the Dirac operator d is thought of as "small" and Eq. 20 is expanded for small d, D ~ d. Feynman rules for 
the Wilson overlap action have been given by Ref. [7], and can be straightforwardly be adapted for any kernel action. 
For overlap actions, it is customary to define the massive overlap operator in terms of the massless one. as 

% 5 ) = (l-^)fl(0)+m, (21) 

This results in an annoying entanglement of the mass with the vertices, complicating a direct computation of the 
running mass. Fortunately, we can compute the multiplicative renormalization factor for the fermion mass indirectly 
as the inverse of the scalar current renormalization factor, and we can evaluate the latter expression at zero quark 
mass. 
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FIG. 1. Three-point and four-point vertices, showing my convention for momentum flows. 

C. Unitary Fat Links 

We will be concerned only with unitary fat links, gauge connections which are themselves elements of the gauge 
group, even though they may be built of sums of products of the original thin links of the simulation. For smooth 
fields the fat links have an expansion VJj(x) = 1 + iaB^x) + . . . and the original thin links have an expansion 
U^(x) = 1 + iaA^(x) + . . .. For computations of 2- and 4-quark operator renormalization/matching constants at one 
loop, only the linear part of the relation between fat and thin links is needed, and it can be parameterized as 

M x ) = h»v{y)A v {x + y) . (22) 

V," 

Quadratic terms in (22), which would only be relevant for tadpole graphs, appear as commutators and therefore do 
not contribute, since tadpoles are symmetric in the two gluons [18,12-14]. In momentum space, the convolution of 
Eq. (22) becomes a form factor 

B^(q) =^h^(q)A I/ (q) . (23) 

V 

The reader could think of fat-link action Feynman rules as being constructed in two levels: First find the vertices for 
actions with ordinary thin links, and then replace the thin link by a unitary fat link. Each quark-gluon vertex gets a 
form factor h^ lv {q) 1 where q is the gluon momentum. If all gluon lines start and end on fermion lines, then, effectively, 
the gluon propagator changes into — ► hn\G\ a h av . 

(Notice that this is a perturbative realization of the statement that fat link fermion connections can be converted 
into thin link fermion connections by redefining the fat link variable as an ordinary thin link variable, but with a more 
complicated pure gauge action.) 
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r 




(c) 

FIG. 2. The three one-loop diagrams (a) "sunset" and (b) "tadpole" fermion self energy and (c) vertex renormalization, 
showing my convention for momentum flows. 

Form factors for APE and HYP-blocked links are given in Appendix B. 

D. One-loop Diagrams 

The "standard method" for carrying out a lattice perturbation theory calculation is to construct the integrand 
I(k,p, m, A, a) analytically as a combination of terms multiplying Dirac matrices, and then to project out the desired 
Dirac structure, producing a single scalar expression I(k,p, m, A, a). Possible divergent terms are subtracted out, and 
the integral is performed using a Monte Carlo routine like VEGAS [20]. The lattice actions I have studied are rather 
complicated (as are most improved actions) and this procedure appears at first sight to be somewhat daunting. The 
difficult part of the calculation is the Dirac reduction, particularly as the action (and its Feynman rules) become 
complicated. One can, of course, do the Dirac algebra using a symbolic manipulation code. However, there is a 
simpler path for the typical lattice practitioner: Take the parts of the programs which one has already written to do 
full-scale numerical simulations of the lattice action, and extract the routines which multiply Dirac matrices times 
spinors. Write the vertices and propagators as explicit 4x4 matrices, string the expressions together, a and let the 
computer do all the Dirac multiplication and projection as if it were doing a standard lattice Monte Carlo calculation. 
All calculations factorize into a set of routines for each action of interest, and a set of integrands (with appropriate 
projection algorithms) for each coefficient. 

This "method" is computationally inefficient, but it is easy to study a wide variety of lattice actions. 

I will parameterize the orientation of fermion and gluon momenta through the vertices as shown in Fig. 1. 

The fermion self energy is parameterized for small fermion momentum p and small mass m as 

S(p,m) = S + i7-pSi+mE 2 - (24) 
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If nonzero, So is (minus) the additive mass renormalization. Si is the wave-function renormalization, needed for 
all external lines in vertex functions. The translation of the quark mass from lattice to MS regularization (Z m ) is 
proportional to the difference Si — S 2 . While S is finite (in lattice units; it is proportional to 1/a and so diverges 
in the continuum limit), Si and S 2 are infared divergent. 

S itself is a sum of two terms, as shown in Fig. 2. The "sunset" graph, Fig. 2(a), uses first order vertices: 

Z a =g 2 f V^(p,p + k,-k)S(p + k)vy\p + k,p,k)G^(k). (25) 

Jk 

The "tadpole" graph, Fig. 2(b), is 

Zb = ~l V$(P,P, k, -k)G^(k). (26) 

Momenta are labeled as shown in the figure. Recall that we usually set p = 0. 

To extract Si and S 2 , I expand the propagators and vertices in a power series in p and m (respectively) and keep 
the leading term. If the free Dirac operator is of the form of Eq. (18) this is straightforward to do. I have performed 
this expansion analytically, since these two quantities are IR divergent, and require subtraction. Below, the quantity 
Si will be used, with where S x = g 2 C F /(16ir 2 )Si. 

The renormalization of currents involves both the vertex graph and Si. I compute the vertex graph simply by 
taking the amplitude 

V T = fv^(p,p + k,-k)}S(p + k)T[S(p + k)V^(p + k,p,k)}G^(k), (27) 

J k 

evaluating it at p = 0, m = (with gluon mass A) and tracing it with the appropriate Dirac projector, before doing 
the integral. For the vector and axial currents, I average over traces in the four cardinal dimensions. 

An alternate parameterization of the vertex allows a connection to matching coefficients of the four fermion operators 
of the effective field theory of electroweak interactions: Write 

V T = K T + K llt ,T^ + K 2lfllv T lvlfl + .... (28) 

Grouping propagator-vertex products, 

V r = f T 1 (p,fc)rT 2 (p,fc)G At ,(fc), (29) 

Jk 

withT 1 ^, fc) = (p,p+k,-k)S(p+k) andT 2 (p,fc) = S(p+k)V^ (p,p+k, -fc). One can find the K's by projecting 
the T's onto elements of the Clifford algebra, T = T + 7^1 + a^T 2 + . . ., with T = l/4Tr T, T x = l/4Tr -y^T, and 

T 2 = 1/8 Tr (7p7, - W„)T. (30) 

Direct computation plus a consideration of lattice symmetries then allows us to extract the separate terms of Eq. 
(28), (here in Feynamn gauge) as 

Ko = f (TqTq - 2T 2 1 T|)G W . (31) 

J k 



and 



= J k \nnc^. (32) 



The Wilson and clover actions have only the listed terms (-Koj K\, K 2 ). The overlap action only has nonzero Kq 
and K 2 terms. The continuum calculation with massless fermions only has nonzero K 2 . Finally, actions which only 
approximate an overlap action could in principle span the Clifford algebra, although the coefficients of the other terms 
would be small if the action were a good approximation. This happens for the non-overlap planar action with HYP 
links. 
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The K\ term makes its presence felt most malignantly in the one loop correction to four fermion operators, where it 
is responsible for "bad" operator mixing into opposite-chirality operators. It poisons lattice calculations of Bk with 
Wilson- type quarks. 

To find the full Z-factor, and q* , we also need the coefficients of J\ and J 2 (in the notation of Eq. 13). These are 
recorded in Table I. These results are certainly not new, but it is useful to collect them. 

If only Kq, K\, and K 2 are nonzero, we can immediately write down relations (which appear many times in the 
literature) between the lattice parts of the matching coefficients of the scalar (S), pseudoscalar (P), vector (V), axial 
vector (A), and tensor (T) currents: The only one we will need below is 

Is - Ip = 2{I V - I A ) = 8KL (33) 

These relations can be used to relate the matching factors for four- fermion operators to those for bilinears. The 
earliest reference I can find for this decomposition is by Martinelli [21] and it has been written down most usefully 
by Gupta, Bhattacharya, and Sharpe [15]. A complication which arises in this case is the prescription used to define 
75 away from four dimensions. The combination of 1/e factors from integrals related to Eq. 2 and e factors from the 
Dirac algebra is different in the bilinear and four-fermion cases, meaning that a particular four-fermion Z— factor into 
a particular continuum convention is a linear combinations of bilinear Z^s plus extra constant terms. 

I found it most straightforward way to find these constants was to do the continuum Dirac algebra using the 
techniques of Rcf. [22] (basically copying the examples of Ref. [23]). In order to extract the momentum scale q* , we 
need to separate the A and B coefficients of Eq. 2. Most lattice calculations do not include penguin graphs, and four 
operators are needed for the most frequently performed four fermion matrix elements, combinations of 

0={q^T iq f)®{q^T 2 qf). (34) 

Special cases are (a) Ti = T 2 = 7^(1 — 75): if color labels a = 5, (5 = 7, O = 0\\ if a = (3, 7 = 5, O = 2 ; and (b)thc 
isospin 3/2 operators for electroweak penguins, 

0? /2 = (s Q 7 M (l - lb)d a )[{u f 3"f^{l +75)^/3) - {dplni 1 +75)^/3)] + (s«7 M (l - 75)u a )(u/37M( 1 + 75)^/3) (35) 

and 

0>l /2 = {saJnO- - 15) dp) [(%7 i u( 1 + l5)u a ) - {dp-fnil + 7b)d a )} + (s q 7m( 1 ~ 75)«/3)(%7 i u( 1 + lh)d a ). (36) 
Ingredients for the non-penguin mixing factors are given in Table II. 

III. SOME EXAMPLES 

I have studied a large variety of lattice actions with thin and fat links. Numerical integrals are checked, when 
possible, by comparison against published results. For currents, I vary the gauge choice £ and check that integrals 
(and integrands) remain gauge invariant. As a general rule using double precision insures that the integrands are 
gauge invariant point by point to a few parts in 10 5 . 

My results for any standard thin link action (Wilson, clover, Wilson overlap, . . . ) are (with one exception) not 
new, and there is no point in republishing old results already in the literature. Results for a standard action with 
HYP links are mostly unpublished, so I will show them in tables. I will show results for the clover and planar actions, 
since the HYP clover action might be an attractive action for simulations. I will also tabulate results for the planar 
action, with thin and HYP-blocked links, and the HYP-blocked planar overlap. 

HYP blocking is characterized by three parameters with "preferred" values a\ = 0.75, ot 2 = 0.6, a 3 = 0.3. While 
HYP blocking is typically presented (and used) at this specific value of these coefficients, it is very useful to show 
how the matching coefficients vary with the degree of fattening of the link. I will do this simply by multiplying the 
standard HYP coefficients by an overall scale factor, and tune the scale factor from zero (corresponding to a thin link) 
to 1 to 1.5. There are many other ways to tune HYP blocking, of course. 

With HYP smearing, the first-order formula Eq. (7) is small because of cancellations of negative and positive 
contributions in the integral, and the second order formula is often needed. The reader will note many "cusps" in the 
q* plots for HYP actions as I switch from first to second order q*'s. In contrast, as a general rule, under increased 
APE-smcaring q* usually falls slowly towards zero, and the first-order formula for q* works well, unless the actual 
matrix element vanishes. 
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A. Additive mass renormalization 



All of the features of fat link perturbation theory can be seen in the additive mass renormalization for non-overlap 
actions. Fig. 3 shows the additive mass renormalization for thin link Wilson and clover fermions, with c$w = 1; an d 
for the planar action with HYP-blocked links. The graphs show So, with the definition 5m — a s (q*)So- We see that 
all these thin link actions have large additive mass renormalization. The addition of the clover term reduces So by 
about half, but it is still big. Smearing the gauge fields has a dramatic effect on So, until the scale factor for HYP 
blocking exceeds unity. At this point the blocking enhances the large gluon momentum region in the integrand rather 
than suppressing it. This effect shows little dependence on the choice of fermion action (planar or clover). 




3.5 



3.0 



2.5 



2.0 




FIG. 3. (a) So parameterizing additive mass renormalization for thin link Wilson (cross), and HYP- link link clover (diamond) 
and planar (octagons) fermions and (b) their momentum scale q* with HYP blocking, with the optimum parameters scaled by 
the shown overall scale factor. 



B. Fermion Bilinears 



Next we turn to results for fermion bilinears, parameterized as Zj = 1 + Zig 2 (q*)CF/(16ir 2 ). Table III shows 
Z-factors for currents for the thin link and HYP-blocked planar action. The values quoted in all tables have an 
uncertainty smaller than ±1 in their rightmost digit. 

Table IV shows Z-factors for the thin link and HYP-blocked naive fermions. These results form a (tiny) subset of 
an extensive calculation of matching factors for staggered fermions by Sharpe and W. Lee [13]. The q* values are 
new. 

Matching factors and q* scales for the local vector and axial vector currents for clover fermions and planar fermions 
as a function of fattening strength are shown in Figs. 4 and 5, and results for the local scalar and pseudoscalar 
currents are show in in Figs. 6 and 7. The qualitative features of fattening are the same for both actions: for thin 
link actions all z's are large in magnitude and the differences between "chiral partners" (such as z v and za) are also 
large. Either HYP action has tiny z's (order unity) with differences an order of magnitude smaller. 
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FIG. 4. (a) zy and za, the coefficients of g 2 Cf /(16tt 2 ) for the local vector (octagons) and axial current renormalization 
constants (crosses) for clover fermions with HYP blocking, (b) Momentum scale q* for the local vector and axial currents. 
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FIG. 5. (a) zy and za and their q*'s, for the HYP-planar action, labelled as in Fig. 4. 
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FIG. 6. (a) zs and zp, the coefficients of g 1 'C 'f / '(16-7T 2 ) for the local scalar and pseudoscalar current renormalization constants 
(evaluated at pia = 1) for clover fermions, labelled by octagons and crosses, respectively, (b) Momentum scale q* . 
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FIG. 7. (a) zs and zp and their q*'s, for HYP-planar action, labelled as in Fig. 6. 

Results for overlap actions parallel those for nonchiral actions: matching factors drop when the clover term is 
included, and drop more when the links are fattened. To illustrate this, I present results for the thin link Wilson 
and clover overlap, both with x = 1-6, and for the HYP planar overlap, in Figs. 8 and 9. A table of results for the 
HYP-planar overlap is given in Table V. 
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FIG. 8. Coefficient of g 2 'C 'f / '(16-7T 2 ) (a) and momentum scale q* (b) for the local vector (and axial vector) currents for overlap 
fermions with the planar action kernel and HYP links, at xq = 1.6. Octagons label the planar action; also shown are the thin 
link Wilson (cross) and clover (diamond) actions. 
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FIG. 9. Coefficient of g 2 CF / (Wn 2 ) (a) and momentum scale q* for the local scalar (and pseudoscalar) currents for overlap 
fermions with the planar action kernel and HYP links, at xo = 1.6. Octagons label the planar action; also shown are the thin 
link Wilson (cross) and clover (diamond) actions. 



C. Four fermion operators 

Some sample results for four fermion operators are shown in Table VI. My results for (ordinary non-overlap) Wilson 
fermions agree with the DRED(EZ) results of Ref. [24] and the NDR results of Ref. [15]. My Wilson overlap results 
satisfy the connection between Z + , Zy,A, and Zp y s of Ref. [15] and Table II. For the special case of radius Xq = 1 they 
agree with a calculation of P. Weisz [25]. They differ by an overall additive factor of 14/3 (the precise value comes 
from Weisz; I have only determined this factor numerically) from the results of Ref. [9]. NDR four fermion matching 
coefficients for the Wilson overlap action for many radii can readily be constructed from any desired operator using 
the tables of bilincars from Ref. [8] and the results of Ref. [15] or Table II (though finding the q* scale will require 
actually doing some integrals). 
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FIG. 10. (a) Coefficient of <? 2 /(167r 2 ) at scale fia = 1 and (b) momentum scale q* for Z+ for matching lattice and NDR 
overlap fermions with the planar action kernel, with xo — 1.6. Octagons label the planar action; also shown are the thin link 
Wilson (cross) and clover (diamond) actions. 

In Fig. 10 I show results for the HYP planar action kernel, as well as for the thin link Wilson and clover overlaps. 
In all cases I set xo = 1.6. There is a large reduction in z + by converting from a Wilson kernel to a clover kernel even 
without fat links (recall that the planar action includes a clover term). Fattening the links further reduces z + . 
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FIG. 11. Z* , the wrong-chirality mixing amplitude, for clover fermions with (a) APE blocking (octagons for a = 0.25 and 
diamonds for a = 0.45) (b) HYP blocking, with the optimum parameters scaled by the shown overall scale factor. 
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Finally, to return to non-chiral actions, we can ask how fattening alters the mixing into different chiral sectors. 
This quantity is parameterized by the coefficient Z* = —SKi. From Eq. 28, Z* — zy — za = 2{zp — z$). As we have 
already seen, fattening pushes all the z's closer to zero, and so their differences also become small. Values of Z* for 
the clover action with APE blocking and HYP blocking are shown in Fig. 11. As expected, either smearing can cut 
Z* by over an order of magnitude. 

The Wilson thin link fermion action value Z* = 9.6 [24] is much greater than even the thin-link clover result. 
Converting to fat links without also turning on the clover term helps, but will not be productive: for example, with 
scale factor unity, the HYP-blocked Wilson action has Z* — 2.29, while including a csw — 1 clover term cuts this 
number to 0.16. 

Results for the operators O7 and 0$ for the Wilson and planar overlaps are shown in Table VII. 

Let's use these results to consider a numerical example of the ambiguities which will afflict our calculation of mixing. 
In all cases we will match the lattice and NDR calculations at a scale qi and run the NDR result to qf = 2 GeV 
using the two-loop evolution equation. We will suppose we are using the Wilson gauge action a coupling at /? = 5.9 
and assume that the inverse lattice spacing is l/a= 1.58 GeV. (These are typical numbers from simulations [28]). A 
standard calculation a la Ref. [5] beginning with the logarithm of the average plaquette yields ctjg-g(qa = 1) = 0.199 
and otjjg(qa = w) = 0.137. We begin with the HYP-planar overlap action and imagine just matching at a scale 
qi = \/a or ir/a, and running up. Then the full matching and running matrix is 



and the result at g, = n/a is 



, 0.979 -0.006 , 

Z(q f , qi )Z( qi )=[ _ QMQ LQ99 ) (37) 



„, n ryf \ ( 0.998 ~0.030\ /qq\ 

Z(q f , qi )Z( qi )=^_ 0M3 1Q12 j (38) 



Next, we compute a "lowest order" q* from Eq. 16. We have q^a = 3.6, q^a — 0.68, and 



7f \7( \ I °- 966 -0.020 , ,„„, 
Z(q f , qi )Z( qi ) = ( _ 00Q9 1 n2 ) (39) 



Notice that in all cases the matrices are nearly diagonal. 

In contrast, the first choice for the Wilson action (just set q*a = 1)) gives 



Z(q f , qi )Z( qi )= [ _^Q33 i q^q ) ' ( 4() ) 



matching at q*a — it and running gives 



Z(q f , qi )Z( gi )= [ } 4 ^ Q ^J, (41, 



while using Eq. (16), with q*a~ (3.0, 2.3), we have 



„, . , . / 1.459 -0.016 \ , . 

Z(q f , qi )Z( qi )=^_ QM2 153Q y (42) 

The off diagonal terms in Z are small and not too dependent on the running scheme, but the diagonal coupling of 
the operators suffers a giant renormalization. One probably should treat perturbation theory results for the Wilson 
overlap cautiously for these factors. 



IV. CONCLUSIONS 



Fattening the gauge connections of a lattice fermion action is a simple way to reduce the size of perturbative matching 
coefficients. HYP blocking is a particularly felicitous choice: it combines large scale smoothing with locality. It is 
clear from the results presented that the qualitative features of fattening do not depend on the specific choice of 
parameters. In situations where full chiral symmetry might not be necessary, the HYP-blocked clover action might 
be an attractive choice for a light quark action. 
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In this work I have only considered lattice actions with scalar and vector couplings and nearest and next-nearest 
couplings. It would be easy to construct the Feynman rules for "hypercubic" actions (such as those of Refs. [19] or 
[26]). Techniques similar to the ones used here could be applied (with only a little more effort) to the more complicated 
approximate fixed point or overlap actions used by several other authors, whose kernels fill out the entire Clifford 
algebra [27]. 

The result for Z + will be used in a lattice calculation of Bk using the HYP-blocked planar overlap [28] . 
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APPENDIX A: FEYNMAN RULES FOR PLANAR-ACTION FERMIONS 

The non-overlap fermion action I am most interested in has scalar and vector couplings to fcrmions offset on 
nearest- neighbor and diagonal-offset sites, and minimal length gauge paths built of unitarized fat links connecting 
them. The nearest-neighbor vector and scalar couplings are labeled p\ and Ai; the diagonal (r = ±/x ± v, v ^ p) 
couplings are p2 and A2. There is also a local scalar coupling Aq. For the Wilson or clover action, p 2 = A 2 = 0; 

pi = Ai = -1/2. For the "planar" action of Ref. [2] Ai = -0.170, pft = —0.177 and diagonal neighbors (r — ±p,±v>, 
(2) (2) 

v 7^ p; A 2 = —0.061, pjx = p v = —0.0538. The constraint A= — 8A1 — 24A 2 enforces masslessness on the free 
spectrum, and —1 = 2p^ + 12 p^ normalizes the action to —ipi-y^dfiip in the naive continuum limit. 
The free fermion action is then 

d(p) = A + 2Ai^cosp AJ +4A 2 ^ cos q^ cos q v + i ^ 7^ sin q^ (2pi + 4p 2 ^ cos q v ) . (AI) 

Nearest-neighbor connections basically contribute rescaled versions of the usual Wilson-action Feynamn rules. The 
diagonal-offset neighbor gauge connections are taken to be an average of the two length-two shortest paths connecting 
the fermions (each of which is a product of unitarized HYP links) . Neglecting the form factor arising from fattening, 

the first order vertex is 

V™ = 2^{(pi + P2 E^(cos Pl + cos P2)v) cos( l£l±pk} 
+ 7,(2*P2 sin((£i^k)(sinp l!y + sinp 2ly )) 

- 2(Ai + A 2 E^(cospi + cosp 2 )„) sin(to^). 

(A2) 

I also include a clover term in V^: 

v^ cl) = -\c swl ^ ^ sinfc - cos y • ( A3 ) 

(My definition of Csw would be unconventional for Wilson fermions if Ai 7^ 1/2, or r 7^ 1 in usual usage.) The planar 
action has Csw — 1-03. 

The expression for the four-point vertex is long, but in all the calculations done here, I only need an expression for 
the vertex at zero fermion momenta and for its first derivative with respect to (equal) fermion momenta, also at zero 
momentum. In that limit, 

V$> = V(-2Ai-12A 2 ) + (l- V)4A 2 sm^sin^) (A4) 

and its derivative at zero external fermion momenta is 

dV^ k k 

-jf- = fVy*(W- 2 Pi - 24 <°2) + (1 - <W 4 P2 sin sin -£) (A5) 



14 



APPENDIX B: EXPLICIT FORMULAS FOR FAT LINKS 



APE blocking: The link after n + 1 smearings is related to the link after n smearings by 

V^ n+1) (x) = Proj su{3) ((l - a)V^\x) 

+ «/6 E^(K W (x)v{r ] (x + u)vS n) (x + HV 

+ V!) n \x - v^V^ix - 0)vS n) (x -v + A)))- (Bl) 

V^ n+1 \x) is projected back onto SU(3) after each step, and v/f \n) = U^n) is the original link variable. The 
momentum-space smearing factor for one level of smearing is 

= /(?)(<W - ^) + ^ , (B2) 

with = |sin(^2ii) and f(q) = 1 — fq 2 . After N smearings, h^ v {q) becomes h^ v {q), which is just h^ v with / 
replaced by f N . 

HYP blocking: The momentum space version of HYP blocking is 
B M = V0?)^0?) + 0(A 2 ), where 

1/ 

p 

n^q) = l + a 2 (l + a 3 )-^(l + 2a 3 )(f-ql-ql) + ^ J] q 2 . 



with ai = 0.75, a-i = 0.6, and <x 3 — 0.3 the favored parameterization of Ref. [3]. 



Ilatt A B 

Z v 4:K 2 + IKx +K -Si 
Z A 4A" 2 - 2K X + K - Si 
Z P I&K2 - 4Jd +K -Si 3 -1/2 
16if 2 + AKi +K0-S1 3 -1/2 



TABLE I. Ingredients for matching coefficients for some local operators. 



Ilatt A B 

Z + 8/Z{K 2 +K - S{) -2 -5/3 

Z_ 80/37^2 - 8/3(7^0 - Si) 4 10/3 

Z 77 4/3(5i<: 2 +2(^o -Si)) -1 7/6 

Z 8S 4/3(32^ 2 +2{K - Si)) 8 -1/3 

^78 12^2 3 -7/2 

Z 87 0-3 



TABLE II. Ingredients for the (non-penguin parts) of matching coemcients for some four-fermion operators. 
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action 


process 


Zi 


* 


thin planar 










So 


3.71 


3.19 




Zv 


-21.5 


2.8 




Za 


-20.0 


2.8 




Z.5 


-26.2 


2.4 




Zp 


-29.3 


2.4 


HYP planar 










So 


0.26 


2.37 




Z v 


-2.27 


1.55 




Z A 


-2.19 


1.54 




Zs 


-0.92 


1.50 




Z P 


-1.11 


1.51 



TABLE III. Table of Z-factors and g*'s for planar and HYP-planar actions, defined so Zi = 1 + Zig 2 Cf / (16tt 2 ). Csw = 1.03. 



action 


process 


Zi 


9* 


thin 










Zv,a 


-14.8 


3.27 




Zp,s 


-39.2 


2.79 


HYP 










Zv,a 


-0.945 


2.36 




Zp,s 


-0.592 


2.05 



TABLE IV. Table of Z-factors and g*'s for thin link and HYP link naive fermion actions. 



process 


z 


q' 


Zv,A 

Zp,s 


-0.489 
0.53 


1.26 
1.96 



TABLE V. Table of Z-factors and q*'s for bilinears, HYP-planar overlap action. 



action 


2+ 


q* 


thin Wilson 


33.4 


2.86 


thin Clover 


14.3 


4.07 


0.45 x 10 APE Clover 


-5.8 


1.05 


0.45 x 10 APE planar 


-5.67 


0.95 


HYP planar 


-3.97 


0.92 



TABLE VI. Table of Z-factors into NDR and q*'s for O+ = Ox + O2 for overlap actions with various kernels. In all cases 
xo = 1.6, Csw = 1 for clover action and 1.03 for planar action. Z + = 1 + a s (q*)/(4n)(z + + 41og(o//)). 
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action 


■-7 

■£77 


^78 


^87 


■-7 

-^88 


thin Wilson Z 


36.29 


-3.51 


-3 


34.85 


thin Wilson L\ 


70.57 


1.56 


-0.003 


75.26 


thin Wilson £2 


17.6 


-4.64 


9.84 


19.32 


HYP planar Z 


-0.65 


-1.98 


-3 


2.41 


HYP planar £1 


1.52 


2.72 


-0.003 


9.69 


HYP planar C 2 


4.94 


-1.65 


9.84 


27.50 



TABLE VII. Table of Z-factors and q*'s for O7 and Og for overlap actions with various kernels. In all cases £0 = 1-6, 
Csw = 1-03 for the planar action. Zij = 5ij + a s (q*) / \Air)Zij . d and Ci are the integrals of Z, weighted by log q 2 and 
log 2 q 2 , respectively. 
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